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ABSTRACT 

We report 31 GHz CARMA observations of IDCS J1426. 5+3508, an infrared-selected galaxy cluster 
at z = 1.75. A Sunyaev-Zel'dovich decrement is detected towards this cluster, indicating a total mass 
of M2oo,m = (4.3 ± 1.1) x 10 14 Mq in agreement with the approximate X-ray mass of ~ 5 x 10 14 M Q . 
IDCS J1426. 5+3508 is by far the most distant cluster yet detected via the Sunyaev-Zel'dovich effect, 
and the most massive z > 1.4 galaxy cluster found to date. Despite the mere ~ 1% probability 
of finding it in the 8.82 deg 2 IRAC Distant Cluster Survey, IDCS J1426. 5+3508 is not completely 
unexpected in ACDM once the area of large, existing surveys is considered. IDCS J1426. 5+3508 
is, however, among the rarest, most extreme clusters ever discovered, and indeed is an evolutionary 
precursor to the most massive known clusters at all redshifts. We discuss how imminent, highly 
sensitive Sunyaev-Zel'dovich experiments will complement infrared techniques for statistical studies 
of the formation of the most massive galaxy clusters in the z > 1.5 Universe, including potential 
precursors to IDCS J1426. 5+3508. 

Subject headings: galaxies: clusters: individual (IDCS J1426. 5+3508) — galaxies: clusters: intra- 
cluster medium — galaxies: evolution — cosmology: observations — cosmology: 
cosmic background radiation 



1. INTRODUCTION 

As the most massive collapsed objects, galaxy clusters 
provide a sensitive probe of the cosmological parameters 
that govern the dynamics of the Universe. Traditional 
cluster cosmology experiments constrain these parame- 
ters by comparing predicted cluster mass fu nctions with 
obser ved cluster counts (for a review, see lAllen et all 

Hp). 

The recent discovery of a handfu l of very massive , 
high redshift (z > 1) galaxy clusters dRosati et al.l 120091 : 
iBrodwin etaD 120101: iFolev et al.l |2011|) has prompted 
several studies of the power of individual massive, rare 
clusters to probe the physics of inflation, in particular 
non-Gaussianity in t he density f l uctua t ions that seed 
structure formation ( Dalai et al.l 120081 : IBrodwin et al.l 
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IParaniape et al.|[2011[ ). These rare, massive clusters were 
all identified from their hot intracluster medium (ICM), 
whether via X-ray emission or from prominent Sunyaev- 
Zel'dovich (SZ) decrements. 

Conversely, the most successful method for discov- 
ering large samples of high redshift galaxy clusters, 
out to at least z = 1.5, is via infrared selection 
(e.g., lEisenhardt et al.l 120081: iMuzzin et all 12008ft . For 
instanc e, the Spitzer/IRAC Shallow Cluster Survey 
(ISCS, lEisenhardt etafl l2008h . which employs a full 
three-dimensional wavelet search technique on a 4.5/im- 
selected, stellar mass-limited galaxy sample, has dis- 
covered 116 candidate galaxy clusters and groups at 
z > 1. More than 20 of these ha ve now been spectro- 
scopically confirmed at 1 < z < 1.5 ([Stanford et al.ll2005t 
l Elston et alJl2006t IBrodwin et al.ll2006t lEisenhardt et al.l 
2008; Brodwin et al.ll201lD . 

The cosmological importance of individual galaxy clus- 
ters depends primarily on their rarity, with the most 
massive clusters at any redshift providing the most lever- 
age. The maximum predicted cluster mass increases with 
time due to the growth of structure. Consequently, clus- 
ters with moderate masses, by present-day standards, 
of a few xlO 14 M Q are cosmologically interesting at 
the highest redshifts (e.g., z > 1.5). O nly one such 
clust er, XDCP J0044.0-2033 at z = 1.58 (jSantos et al.l 
I2011D . has been found to date, although at such a high 
redshift its Lx _ based mass is quite uncertain. The 
most sensitive, large angle, ICM-based cluster survey to 
date, the 2500 deg 2 S outh Pole Telescope Survey (SPT; 
ICarlstrom et al.ll2011l) , would only be able to detect such 
high redshift clusters if they were more massive than 
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«3x 10 14 Mq ([Vanderlinde et alJl2010HAndersson et al.l 
1201 ID . A population of such high-redshift clusters, if 
massive enough, could provide a powerful opportunity 
to study the physics of the early Universe and, in par- 
ticular, to either confirm or falsify the predictions of the 
ACDM paradigm. They would also be invaluable probes 
of the early formation of massive galaxies in the richest 
environments. 

To extend the ISCS infrared-selected cluster sample to 
the z > 1.5 regime, we carr ied out the Svitze r Deep, 
Wide-Field Survey fSDWFS. lAshbv et al.l [2001) . which 
quadrupled the Spitzer/IRAC exposure time, leading to 
a survey more than a fa ctor of 2 deeper than th e original 
IRAC Shallow Survey ([Eisenhardt et all 120041) . A new 
cluster search, the Spitzer/IRAC Distant Cluster Survey 
(IDCS), to be described in detail in forthcoming papers, 
was carried out using new photometric redshifts calcu- 
lated for the larger, deeper SDW FS sample. Cluster 
IDCS J1426.5+3508 at z = 1.75 (|Stanford et al.l [2012L 
hereafter S12) is the first of the new high-redshift IDCS 
clusters to be spectroscopically confirmed. 

In this paper we report the robust detection of the 
SZ decrement associated with IDCS J1426. 5+3508 in 
31 GHz interferometric data taken with the Sunyaev- 
Zel'dovich Array (SZA), a subarray of the Com- 
bined Array for Research in Millimeter-wave Astron- 
omy (CARMA). This is by far the most distant clus- 
ter yet detected in the SZ effect, a consequence of the 
fact that IDCS J1426. 5+3508 is the most massive clus- 
ter yet discovered at z > 1.4. In ij2]we briefly describe 
the discovery of cluster IDCS J1426. 5+3508. In g3] we 
present the CARMA observations and reductions, and 
describe the measurement of the Comptonization and 
mass of IDCS J1426. 5+3508. We also discuss various 
systematic uncertainties present in our analysis. In [jH 
we calculate the likelihood of having found this clus- 
ter, given our survey selection function, in a standard 
ACDM cosmology. In SJH we place IDCS J1426.5+3508 
in an evolutionary context, finding that it is a pre- 
cursor to the most massive known clusters at all red- 
shifts. In [J6] we discuss our results in light of current 
and upcoming SZ surveys. We present our conclusions 
in SJ7] In this paper, we adopt the WM AP7 cosmology 
of jnA ,fl M ,h) = (0.728,0.272,0.704) of lKomatsu etafl 
(|2011[) . Except where otherwise specified (in ij2] and 
§3.2|) . we use a cluster mass (M2oo,m) defined as the mass 
enclosed within a spherical region of mean overdensity 
200 x Pmean, where Pmcan is the mean matter density on 
large scales at the redshift of the cluster. 

2. CLUSTER IDCS J1426. 5+3508 

IDCS J1426.5+3508, first reported in S12, was discov- 
ered in the IDCS as a striking three-dimensional over- 
density of massive galaxies with photometric redshifts 
at Zphot ~ 1-8 (Figure [TJ upper right panel). Follow- 
up spectroscopy, with HST/WFC3 in the infrared and 
Keck/LRIS in the optical, identified seven robust spec- 
troscopic members at z = 1.75 and ten additional likely 
members with less certain redshifts. The HST data also 
revealed the presence of a strong gravitational arc, the 
analysis and i mplications of which ar e discussed in a com- 
panion paper (|Gonzalez et al.|[20l2l ). 

IDCS J 1426 5+3508 is also dete cted in shallow Chan- 
dra data from lMurrav et ahl (|2005| ). with an X-ray lumi- 



nosity of io.5-2koV = (5-5 ± 1.2) x 10 44 erg /s . Usin g the 
Mr, np-Lx scaling rela tion of iVikhlinin et al.l ([20091 ) and 
the lDuffv et al.l ([2008D mass-concentration relation, this 
corresponds to a mass of M 200 c =2 5.6 x 10 14 M (S12J3. 
The HST/ ACS and WFC3 color image, along with the 
X-ray detection (yellow contours) are shown in the lower 
panel of Figure [j 

While IDCS J1426.5+3508 appears from the X-ray 
data to be impressively massive, particularly at such an 
extreme redshift, the X-ray luminosity is based on only 
54 cluster X-ray photons. Further, it relies on a high- 
scatter mass proxy, the Mqqq jC -Lx relation, which is only 
calibrated at low (z < 1) redshift. The uncertainty of this 
X-ray mass estimate is at least ~ 40%. 

3. DATA 

3.1. Sunyaev Zel'dovich Observations 

IDCS J1426. 5+3508 was observed with the SZA, an 8- 
element radio interferometer optimized for measurements 
of the SZ effect. During the time in which the data were 
obtained, the array was sited at the C edar Flat location 
described in iCulverhouse et all (|2010| ) , and was in the 
"SL" configuration. In this configuration, six elements 
comprise a compact array sensitive to arcminute-scale SZ 
signals, and two outrigger elements provide simultaneous 
discrimination for compact radio source emission. The 
compact array and outrigger baselines provide baseline 
lengths of 0.35-1.3 kA and 2-7.5 kA, respectively. 

The cluster was observed for four sidereal passes on 
nights ranging over UT 2011 June 9-29, for a total 
of 6.24 hours of unflagged, on-source integration time. 
Data were obtained in an 8 GHz passband centered at 
31 GHz. The resulting noise level achieved by the short 
(SZ-sensitive) baselines was 0.29 mJy/beam, or 22 /iKrj 
with a synthesized beam of 118" x 143". For the long 
baselines, the noise level achieved was 0.29 mJy/beam 
in a 19'.'9 x 14'.'4 beam. The data were reduce d us- 
ing the Miriad software package ([Sault et al.lll995|) . and 
the absolute calibration was determined us ing periodic 
measu rements of Mars calibrated against the lRudv et al.l 
(1987) model. A systematic uncertainty of 5% has been 
applied to the SZ measurements due to the uncertainty 
in the absolute calibration. The SZA map is shown in 
the upper left panel of Figure [TJ 

3.2. Comptonization and Mass of IDCS J1426. 5+3508 

To determine the properties of IDCS J1426. 5+3508 
from the SZA data, we make use of the fact that the 
integrated Compton y-parameter scales with total mass. 
Assuming a spherically symmetric pressure distribution 
within the cluster, the integrated y-parameter is given 
by 

YD\ = ^-J P(r)dV, (1) 

where Da is the angular diameter distance, o~t is the 
Thomson cross-section of an electron, m e is the electron 
mass, c is the speed of light, P(r) is the pressure profile, 
and the integral is o ver a spherical volum e centered on 
the cluster (see, e.g., IMarrone et al]|201lD . The pressure 

13 In S12 the mass is measured with respect to a critical over- 
density, rather than mean overdensity we use in this paper. 
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Figure 1. Upper left: 12.8' X 12.8' SZA mm-wave map showing the SZ decrement of IDCS J1426. 5+3508. The flux scale in mjy/beam is 
shown to the right of the map, and the FWHM of the CLEAN beam is shown in the lower left corner. Upper right: A 5' X 5' By/R[A.S\ 
pseudo-color image of IDCS J1426. 5+3508, with the SZ contours overplotted. The optical and infrared data are from the NOAO Deep, 
Wide-Field Survey (NDWFS: Dannuzi and Dcv 1999) and SDWFS, respectively. Lower panel: HST F814W+F160W pseudo-color zoom-in 
of IDCS J1426. 5+3508. The SZ contours (green) from the present work and X-ray contours (yellow) from S12, corresponding to 0.131, 
0.077 and 0.044 counts per 2" X 2" pixel in 8.3 ks in the 0.5-2 KeV band, are overplotted. Two sources are indicated by hash marks. The 
northern source is an X-ray luminous QSO that is also a cluster member; the southern source is a radio-bright AGN, the emission from 
which was removed from the SZA map. Both are discussed further in the text. In all panels North is up, East is left, and a scale bar is 
given. 
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profile P(r) is determined by fitting a m odel profile to 
the SZ A data. We adopt the generalized iNavarro et al.l 
(|1997D model proposed bv lNagai et~aT1 ([2007T ): 



P(r) = 



Pa 



xT(l + a; Q )^-T)/ Q ' 



(2) 



where x — r/r s . We fix the power law expo nents (ct, /3, 7) 
to th e "universal" values reported b y lArnaud et all 
(2010f), and leave Pq and r s free to varjQ- 

One important source of systematic errors in SZ cluster 
measurements is radio point sources. These synchrotron 
sources are often variable, are known to be correlated 
with clusters, and c an fill in the SZ decrem ent at low fre- 
quencies (see, e.g.. iCarlstrom et al.l [2002h . Wide angu- 
lar scale surveys such as the SPT survey generally mask 
regions where bright point sources are detected, which 
changes the survey selection function but can be ac - 
counted for using simulations (jVanderlinde et all 120101 ) . 
However, for targeted SZ follow-up of sources selected in 
other bands, point source fluxes must be measured and 
removed in order to accurately estimate the cluster mass. 
Due to the temporal variability of these sources, the point 
source fluxes ideally should be determined concurrently 
with the SZ measurement. The SZA uses its longer 
outrigger baselines to perform such simultaneous point 
source measurements. For spatially compact sources, the 
emission measured on the outrigger baselines (smaller 
angular scales) can be assumed to be representative of 
the emission measured on the short SZ-sensitive base- 
lines (larger angular scales) , and can thereby be rem oved 
from the SZ signal (see, e.g.. iMuchovei et al.ll2007l ). 

Two such point sources were discussed in S12 and are 
indicated by hash marks in Figure Q] The northern 
source is an X-ray bright QSO that is also a spectro- 
scopically confirmed member of IDCS J1426. 5+3508. As 
described in S12, the contribution to the cluster's X-ray 
flux from this QSO had to be removed to obtain a reliable 
cluster mass. This QSO is not detected in the radio, how- 
ever, with a formal 31 GHz flux density of —0.04 ± 0.31 
mjy, and therefore it has no impact on the measured 
cluster mass. A second point source, bright enough in the 
radio to fill in the SZ decrement of IDCS J1426. 5+3508, 
was identified near the cluster center at 14 h 26 m 32!l, 
35°08'15'.'2, indicated by the southern hash mark in Fig- 
ure [TJ This source corresponds to a ~ 95 mjy source 
identified at 1.4 GHz in the NV SS (|Condon et al.|[l99ll 
and FIRST (|Becker et al.l 119951 ) catalogs. The lower- 
frequency radio survey data indicate that this source is 
unresolved at the angular scales to which the SZA is 
sensitive, allowing it to be modeled as a point source 
and constrained using the long baseline data. The flux 
density of this source is left as a free parameter in the 
cluster model fit, and is found to be 5.3 ± 0.3 mjy. The 
resulting spectral index of —0.93 + 0.02 is consistent with 
synchrotron radiation. 

A Markov Chain Monte Carlo (MCMC) fit is used to 
jointly constrain Pq, r Sl and the flux densities of the two 
point sources in our fiducial cosmology. The results of 
this fit are used to compute the integrated y-parameter 



14 lArn aud ct al. (2010) suggest other power law exponents ap- 
propr iate for clusters with a known morphology, but Marr one et al.l 
(2011) find that these alternative choices make a < 5% difference 
in the estimate of Ysoo- 



and the total cluster mass. By convention, both the total 
cluster mass and the integrated y-parameter are reported 
within a radius corresponding to a fixed overdensity 
A = 500 relative to the critical density of the Universe 
at the redshift of the cluster. The integrated quantities 
are determined by finding the overdensity radius r§oQ 
(equivalent to the mass, since Ma = (47rr A /3)Ap c (z)) 
and integrated ^-parameter Y500 that a re mutually con- 
sistent with the I Andersson et all ([201 ID M 500 -Y 50 o scal- 
ing relation. We choose this relation since the clus- 
ters with which it is measured have the highest mean 
redshift. The self-consistent values of rsoo and Y500 
are determined iteratively at each step in the MCMC, 
and the scaling relation parameters are sampled at each 
step using their reported means and variances. Using 
this method, we find that M 50 o,c = (2.6 ± 0.7) x 10 1 * 
Mq , where the error i s statistical in nature. Assuming 
the IDuffv et al.l (|2008f ) concentration relation, this cor- 
responds to M 200iC = (4.1 ± 1.1) x 10 14 M Q . This mass 
agrees well with the X-ray mass of Af 2 oo,c — 5.6 x 10 14 
Mq, for which the error is at least 40%. 

A significant uncertainty, inherent to both SZ and X- 
ray measurements of the masses of high redshift clus- 
ters, is due to the fact that no published mass-observable 
scaling relation has been estimated using a sample con- 
taining cluster s at redshifts significan tly greater than 
unity. While [Andersson ct al. (2011) found no com- 
pelling evidence for redshift evolution in the M500-I500 
scaling relation out to z ~ 1, their sample consisted 
of lower-redshift and higher-mass objects whose scaling 
properties may be systematically different from those of 
IDCS J1426. 5+3508. 

To enable the SZ mass of IDCS J1426. 5+3508 to be 
easily calculated with future scaling relations, we also re- 
port the spherically-averaged dimensionlcss Comptoniza- 
tion, y sp h,5oo = (7-9 ± 3.2) x 10~ 12 . This corresponds to 
a gas mass at r5oo of A/ gas = (2.9 ±1.2) x 10 13 M© for a 
constant temperature of 5 keV. This temperature is con- 
sistent with the X-ray measurement (S12); for a cooler 4 
keV cluster the inferred gas mass would increase by 25%. 

To assess the SZ detection significance, we fit a clus- 
ter model to the visibility data and compare to a null 
model. Our model consists of an elliptical Gaussian SZ 
decrement located at the X-ray centroid along with the 
nearby radio point source, and our null model consists 
only of the point source. We find a A\ 2 of 37.65 between 
the model and null model fits. Since the (fixed-position) 
elliptical Gaussian adds four degrees of freedom, this A% 2 
corresponds to a detection significance of 5.3 a, or a false 
detection probability of 1.32 x 10 -7 . Because the SZ sig- 
nal was sought only at the position of the X-ray emis- 
sion, this detection is very secure. In contrast, a false 
5.3 sigma detection is far more plausible in a wide-field 
survey like SPT because there are of order 10 6 indepen- 
dent pixels in which such a deviation can arise by chance. 
We therefore consider the pointed, follow-up detection of 
IDCS J1426. 5+3508 to be robust. 

If the cluster position is not fixed to the X-ray cen- 
troid, the best-fit centroid of the SZ signal is 14 h 26 m 34!0, 
35°08 / 02'/6, which is slightly offset (« 32" and 25") from 
the X-ray and BCG positions, respectively. As the uncer- 
tainty in the SZ centroid is « 35", this offset is not sta- 
tistically significant. Absent the X-ray positional prior, 
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the significance of the SZ detection is almost identical 
(5.2 a). 

4. PROBABILITY OF EXISTENCE IN A ACDM UNIVERSE 

While not as massive as the z ~ 1 SPT clusters, 
IDCS J1426. 5+3508, by virtue of its exceptional redshift 
and substantial mass, is arguably one of the more ex- 
treme systems discovered to date. It is therefore inter- 
esting to calculate the probability of discovering it in the 
IDCS in a standard ACDM Universe. 

We use the ex c lusion curve formalism presented in 
iMortonson et alJ (|201lL hereafter Mil), which tests 
whether a single cluster is rare enough to falsify ACDM 
and quintessence at a given significance level, account- 
ing for both sample and parameter variance. The Mil 
formalism uses cluster masses measured within a sphere 
defined by an overdensity 200 times the mean, rather 
than critical, density of the Univ erse. In this convention, 
adopting the iDuffv et alJ (|2008[ ) mass-concentration re- 
lation, IDCS J 1426. 5+3508 has a mass of M 20 o,m = 
(4.3 + 1.1) x 10 14 M Q . 

The Mil formalism requires that observed cluster 
masses be corrected for Eddington bias caused by the 
steep slope of the cluster mass function at this mass 
scale and redshift. We follow the prescription given in 
Appendix A of Mil: 

In M = In M obs + ^7°f n u-, ( 3 ) 

where M is the expected value for the true mass given the 
observed mass, 7 is the local slope of the mass function, 
and a is the uncertainty associated with the observations. 
Here we use the fitting function for 7 provided in Mil, 
which yields 7 = —5.8. We thus derive that the expected 
value for this bias-corrected mass is M2oo, m = (3.6 ± 
0.9) x 10 14 Mq. 

In Figure [5] we plot the appropriately bias-corrected 
mass for IDCS J1426. 5+3508 along with 95% con- 
fidence exclusion curves using the Mil prescription. 
For comparison we include rare, massive clusters from 
the S PT survey ([Williamson et all 120111 iBrodwin etafl 
I201C0 and other massive z > 1 clusters from the ISCS 
for which weak lensing or X-ray ma sses are available 
(jBrodwin et all l20ll Uee et alJ l20ll . All masses are 
deboosted in accordance with the Mil prescription. 

The lowest curve in the figure corresponds to the 95% 
exclusion curve for clusters within the 8.82 deg 2 IDCS 
area. The central value of the deboosted SZ mass is 
above the curve, formally falsifying ACDM at the 95% 
level. We note, however, that the statistical error bar 
overlaps the exclusion curve. Furthermore, uncertainties 
in the scaling relation could easily lower the mass below 
the exclusion curve. 

Rather than asking whether such a cluster should have 
been detected in the IDCS, perhaps a fairer question 
is to ask whether such a cluster should have been de- 
tected in the full area of sky covered by all cluster 
surveys capable of detecting such a cluster. Numer- 
ous X-ray surveys of varying areas and depths could 
have detected IDCS J1426. 5+3508, but for simplicity 
we consider the SPT survey, which has the sensitiv- 
ity to det ect clusters like IPCS J 1426. 5+3508, even at 
z = 1.75 (jVanderlinde et al.|[20Toh . The first 178 deg 2 
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Figure 2. Mil-style plot showing the mass and redshift of 
IDCS J1426. 5 +3508 (large red circle), along with other z > 1 
ISCS clusters IBrodwin et al.H2011k Uee et alJl20Tll. small red cir- 
cles) and SPT cluste rs ijWilliamsoii ct al. 2011, small blue circles; 
IBrodwin et al.l 120101 . small cyan circle) . The solid red curve is 
the 95% exclusion curve for the IDCS area. The cyan, blue and 
black curves are the exclusion curv es for the currently pub lished 
full depth SPT survey (178 deg 2 , [Vandcrlindc ct al. 2010]), the 
2500 deg 2 SPT preview survey of Williamson ct al. (2011) and 
the full sky, respectively. The square symbol represents clus- 
ter SPT-CL J0102-491 5, first reported as ACT-CL J0102-4915 
(IMenant cau ct al. 2010). We plot t his cluster at the spectro scopic 
redshift of z = 0.870 reported in IMenanteau et alJ 1120121 ). All 
clusters are color-coded to the appropriate exclusion curve. 

of the SPT survey does not contain a single cluster 
comparable to IDCS J14 26. 5+3508 in mass and redshift 
(|Vanderlinde et all 120101 ) . which suggests that the true 
abundance of such clusters is far lower than our dis- 
covery would indicate. Neither the Atacama Cosmol- 
ogy Telescope (ACT) project nor the Planck project sur- 
yeys are sensitive enough to detect IDCS J1426.5+350 8 
([Marriage et al. 2011; Planck Collaboration et al.ll2011| ). 
If we consider IDCS J1426. 5+3508 to have been drawn 
from an area of at least ~ 178 deg 2 (shown as the cyan 
curve in Fig. [2]) rather than the IDCS volume, then the 
cluster is no longer in contradiction with ACDM; we were 
simply somewhat lucky to detect it in the comparatively 
tiny (8.82 deg 2 ) IDCS area. 

5. EVOLUTION TO THE PRESENT DAY 

Although the mere existence of IDCS J1426.5+3508 
may not have significant cosmological implications, it 
is nevertheless among the rarest, most extreme clusters 
ever discovered. As such, it is interesting to understand 
its nature in the context of the growth of the largest 
structures. 

Using the iTinker et al] (2008) mass function, we 
first identify the space density of clusters at z = 
1.75 that have masses greater to or equal to that of 
IDCS J 1426. 5+3508. At each redshift we then asso- 
ciate the clusters with that space density with its de- 
scendents. In Figure [3] we plot the mass evolution of 
IDCS J1426.5+3508 (large red circle) from z = 1.75 to 
the present day (thick black line). For comparison we 
also plot a representative selection of the most extreme 
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clusters found to date at all redshifts. This abundance - 
matching calculation is designed to predict the mean ex- 
pected growth of IDCS J1426.5+3508. We note that 
there are intrinsic variations in accretio n histories due 
to th e stochastic nature of mergers (e.g., Wcchslcr et al. 
2002), but consider a full treatment beyond the scope of 
the present work. 

As demonstrated in this figure, IDCS J1426.5+3508 is 
much more than just another massive, very high redshift 
cluster. It is an evolutionary precursor to the most mas- 
sive known clusters at all redshifts. Indeed, its descen- 
dent population consists of the most massive virialized 
objects ever discovered. 

There are no other confirmed clusters at z > 1.75 
with masses greater than 10 14 M®. Between 1.5 < z < 
1.75 only one cluster, XDCP J0044 .0-2033 at z = 1.58 
([Santos et al.H2011tlFassbender et al.H201lD . is compara- 
ble to IDCS J1426.5+3508, though it is somewhat less 
massive and below the growth trend. 

The most massive z > 1 X-ra y selected c l uster , 
XMMU J2235.3-2 557 at z = 1.39 (jMullis et all [20051: 
IRosati et al.l 12009). falls along the evolutionary path of 
IDCS J1426. 5+3508, as does the most massive z > 1 
cluster in the SPT s urvey, SPT-CL 2106-584 4 at z = 1.13 
(|Folev et all f20Tll IWilliamson eTaTl l20TTI) . Similarly, 
the only z ~ 1 Planck cluster reported to date, PLCK 
G266.6-27.3 at z = 0.94 is also consistent. 

SPT-CL J0102-4915, first rep orted as ACT-CL J0102- 
4915 from the AC T survey (jMenanteau et all 120101 : 
iMarriage et al.ll20lTh . is among the most significant de- 
tections overall in both of these SZ surveys. Masses 
of this cluster from the two SZ surveys are plotted 
as the squar e symbols at the spect roscopic redshift of 
z = 0.870 (jMenanteau et al.l |20H . The lower SPT 
mass is measured from the SZ decrement, whereas the 
higher ACT mass (slightly offset in redshift for clar- 
ity) is combined from several methods. The masses 
are consistent within the errors. IDCS J1426. 5+3508 
is consistent with gr owing into a cluster like this as 
well. RX .11347-1145 (jAllen et aLllMfil ILu et al.l[20Tfl 
is the most luminous X-ray cluster in the sky, though 
it is not quite massive enough to be considered a de- 
scended of IDCS J1426. 5+3508. In the local Universe 
the Coma cluster represents the kind of extreme cluster 
I PCS J1426.5+3508 wi ll evolv e into. 

Holz and Perlmutterl 1)20101 ) recently predicted the 
properties of the single most massive cluster that can 
form in a ACDM Universe with Gaussian initial con- 
ditions and where the growth of structure follows the 
predictions of general relativity. They predict this clus- 
ter should be found at z = 0.22 and have a mass of 
M 20 o,m = 3.8 x 10 15 M Q . Their 68% contour in pre- 
dicted mass and redshift is shown as a light blue contour 
and labeled HP10 in the figure. IDCS J1426.5+3508 is 
completely consistent with being an evolutionary precur- 
sor to the most massive halo predicted to exist in the 
observable Universe. 

6. DISCUSSION 

Given its extreme mass it is not surprising that 
IDCS J1426.5+3508 was detected in shallow (5 ks) Chan- 
dra images (S12), despite its high redshift. The ro- 
bust, high-significance SZA detection reported herein 
confirms that IDCS J1426. 5+3508 is one of the most 



massive collapsed structures evolving in the Universe. 
Although it was discovered in a small (~ 9 deg 2 ) 
infrared-selected cluster survey, the odds of finding such 
a rare, massive cluster were quite l ow, 0.8% (0.2%) us- 
ing the iHolz and Perlmutterl (|2010l ) formalism with the 
deboosted (measured) mass. Finding the precursors of 
IDCS J1426. 5+3508 at even earlier times will require 
much larger surveys. 

In this section we examine why the current genera- 
tion of SZ surveys, despite covering areas ranging from 
10 2 -10 4 deg 2 , have not yet seen a cluster as extreme as 
IDCS J1426. 5+3508 at z > 1.5. We also discuss the po- 
tential of the next generation of SZ surveys to discover 
its evolutionary precursors. 

6.1. Current SZ surveys 

The detectability of a cluster like IDCS J1426.5+3508 
in an SZ cluster survey depends on three factors. The 
survey must have sufficient mass sensitivity at high red- 
shift, it must cover sufficient area so that the expectation 
value of the number of clusters is greater than unity, and 
the clusters in question must not contain bright radio 
point sources nor have any along the line of sight. 

In terms of sensitivity, only the SPT survey is cur- 
rently capable of detecting IPCS J 1426 5+3508 at this 
high redshift. iVanderlinde et all ((2010) report a for- 
mal 50% completeness of M 2 oo,m ~ 3.5 x 10 14 M Q at 
z = 1.5, and this improves to > 90% at z = 1.75. That 
mass limit is below both the measured and deboosted 
mass of IPCS J1426.5+3508. So the SPT should nom- 
inally have detected such a cluster were it present at 
z = 1.75 in the single-ban d analysis of their first 178 deg 2 
(|Vanderlinde et afll2010D . However, since the SPT com- 
pleteness lim it at high redshift is un certain at the 30% 
level in mass fVandcrlin de et aT1l2010f) . the non-detection 
of a single high-redshift cluster like IPCS J1426.5+3508 
is unremarkable. 

A more interesting question, related to the second cri- 
terion and to the rarity analysis of SjH is that °f now 
many such clusters are actually expected in the cur- 
rent S PT area? According to the IHolz and Perlmutterl 
(2010) formulation, we expect only 0.2 such clusters 
in 178 deg 2 , and therefore none should have been 
fou nd. In the complete, f ull-dep th 2500 deg 2 SPT sur- 
vey, IHolz and Perlmutterl (|2010l ) predict the presence of 
2.4 such clusters, though shot noise will play a role in 
whether any are found. 

Finally, as discussed in g3JJ IPCS J1426.5+3508 has 
a radio-bright AGN at z = 1.535 along the line-of-sight, 
indicated by the southern hash marks in Figure [TJ At 
a redshift separation of Az > 0.2, it is not physically 
associated with the cluster. Its radio emission was con- 
strained using the long SZA baselines, which are insen- 
sitive to the arcminute angular scale SZ decrement. If 
the measured 1.4 to 31 GHz spectral index of this source 
(-0.93) is valid to 150 GHz, the AGN flux would be 1.2 
mjy at this frequency, well below the 5 mJy threshold 
for masking individual point sources in the SPT survey. 
However, it would partially fill in the ~ 14 mJy SZ decre- 
ment at 150 GHz, reducing the SPT detection signifi- 
cance and biasing the mass estimate low by ~ 10 — 15%. 

It is possible that the AGN content of clusters at the 
highest redshifts may be substantially larger than at in- 
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Figure 3. Predicted mass growth of IDCS J1426. 5+3508 vs. redshift based on abundance matching. IDCS J1426. 5+3508 is the large 
red circle and the predicted growth in its mass is shown as a thick black line. The 1 a errors, stemming from the mass measurement 
errors, are shown as the shaded region. An assortment of the rarest, most massive clusters found at all redshift is sh own for comparison 
and d iscussed in the text. SPT-CL J0102-4915 was first reported as ACT-CL J0102-4915 bv lMenanteau et al.l H2010T ) and IMarriage et al.l 
lt20TTT) . The independent mass measurements from both surveys are plotted as the square symbols for this cluster at the spectroscopic 
redshift (z = 0.8701) reported in Mcnantcau et al. (2012). IDCS J1426. 5+3508 is consistent with being a member of the most extreme 
population of virialized structures in the Universe. 



terme diate (z < 1) redshifts. Indeed, Galamct z et al.l 
(2009) find evidence for a strong increase in the inci- 
dence of AGN in clusters with redshift from 0.2 < z < 
1.4, particularly for X-ray and IR-selected AGN. The 
evidence for rapid number density evolution in radio- 
selected AGN is not a s conclusive. In some surveys (e.g., 
iGalametz et al.ll2~009D it increases as quickly in cluste rs as 
in the field, whereas others (e.g.. iGralla et al.ll2011[ ) find 
it is roughly constant. One major advantage of CARMA 
for targeted SZ observations of clusters is the ability to 



simultaneously identify and remove contamination from 
such (generally variable) sources, as in this work, and 
thus recover accurate cluster masses. 

6.2. Next- Generation SZ surveys 

Several next-generation SZ experiments designed 
for CMB polarizat ion measurem ents, such as SPTpo l 
(jBleem et al.l [20121) and ACTpol (jNiemack et al.ll2010l) . 
are currently being deployed. These experiments will 
possess typical SZ mass sensitivities a factor of sev- 



Brodwin et al. 



eral lower than present surveys. Here we calculate 
the expected yield of high-redshift clusters similar to 
IDCS J1426.5+3508 for a canonical 500 deg 2 next 
generation survey with constant M2oo,m completeness 
limits of 1.2 x 10 14 Mq and 1.5 x 10 14 M Q at the 50% 
and 90% levels, respectively. 

The same methodology that explains the lack of a 
cluster like IDCS J1426. 5+3508 in the current SPT 
survey area predicts that a significant number of very 
high redshift clusters, somewhat less massive than 
IDCS J 1426. 5+3508, will be found with the polarization 
experiments. Figure [4] shows the expected c umulative 
yield, based on the lHolz and Perlmutterl (2010) formula- 
tion, at the two mass limits given above, accounting for 
the expected completenesses. 

The combination of area and depth of the next genera- 
tion SZ surveys will allow them to identify a large statis- 
tical sample of massive, hot clusters at 1.5 < z < 2.0. 
As such they provide a natural extension of the in- 
frared surveys, which have characterized galaxy clusters 
to z = 1.5. A 100 deg 2 region that will be covered by 
both SPTpol and ACTpol will also contain Spitzer in- 
frared data, from an ongoing Spitzer Exploration Class 
program (PID 80096, PI Stanford). This will allow opti- 
mal high-redshift cluster detection and study at both IR 
and mm wavelengths, the only two cluster probes that 
are nearly redshift-independent at 1 < z < 2. 

Despite the wealth of new, high redshift clusters to 
be found by the next-generation SZ surveys, there are 
only even odds (p ~ 0.5) of finding az = 2 precursor to 
IDCS J1426. 5+3508. Unless, of course, we are fortunate 
once again. 
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Figure 4. Predicted cumulative counts of massive 1.5 < z < 



2 clusters at the 50% (M2OO 1 



1.20 X 10 14 Mq) and 90% 



(M200 rrl = 1.5 X 10 14 Mq) completeness limits of a representa- 
tive 500 deg 2 next— generation SZ survey. The curves account for 
losses due to incompleteness. The arrow indicates the redshift of 
IDCS J1426.5+3508. 

7. CONCLUSIONS 

We have presented a robust 5.3 a SZ detection of the 
infrared-selected cluster IDCS J1426.5+3508 at z = 1.75, 



by far the most distant cluster with an SZ detection to 
date. Using the lAndersson et ail ()201lD M500-Y500 scal- 
ing relation, we find M 2 oo,m = (4.3 ± 1.1) x 10 14 Mq, in 
agreement with the X-ray mass reported in S12. 

The chance of finding IDCS J1426. 5+3508 in the ~ 
9 deg 2 IDCS was < 1%, which, taken at face value, 
marginally rules out ACDM. Within the much larger to- 
tal area probed by other SZ and X-ray surveys, its exis- 
tence is not unexpected. 

In a ACDM cosmology, with the growth of structure 
as predicted by general relativity, IDCS J1426. 5+3508 
will grow to a mass of A^oo.m ~ 4.7 x 10 15 Mq by the 
present day. It therefore shares an evolutionary path 
with the most massive known clusters at every redshift. 
Along with several of the rarest, high-redshift clusters 
currently known, IDCS J1426. 5+3508 is consistent with 
evolving into the single most massive cluster expected to 
exist in the observable Universe. 
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